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ABSTRACT 


This  Technical  Report  supersedes  TR  No.  7  entitled:  INVERSE 
COMPUTATION  FOR  LONG  LINES:  A  NON-ITERATIVE  METHOD  BASED  ON  THE 
TRUE  GEODESIC,  which  is  out  of  print.  It  contains  the  material  of 
the  original  publication  of  1950,  and  in  addition,  formulas  per¬ 
taining  to  long  lines,  derived  through  the  years  at  AMS. 

The  solutions  of  the  Direct  and  Inverse  Geodetic  Problem  are 
presented  in  forms  which  are  adaptable  to  desk  calculator  and  to 
electronic  computer. 

The  maximum  errors  in  the  solutions  due  to  the  omission  of 
higher  order  terms  have  been  determined  and  are  presented  in  tables 
in  the  Appendix.  These  tables  will  enable  the  user  of  the  solutions 
to  decide  whether  the  accuracy  requirements  can  be  obtained  with  or 
without  the  higher  order  terms.  These  higher  order  terras  have  been 
derived  and  are  presented  herein. 


DIRECT  AND  INVERSE  SOLUTIONS  OF  GEODESICS 


SECTION  I .  GENERAL 

1.  Purpose  and  Scope.  The  purpose  of  thiB  report  is  to  present  in 
a  single  publication  the  various  forms  of  the  Direct  aDd  Inverse 
Solutions  of  Geodesics  which  have  been  solved  by  the  Army  Map  Service. 
This  report  supersedes  AMS  Technical  Report  No.  7. 

SECTION  II.  INTRODUCTION 

In  Section  III  of  this  report  a  procedure  is  given  for  a  rigor¬ 
ous  and  rapid  non-iterative  inverse  solution  of  very  long  geodesics. 
This  procedure,  which  is  in  a  convenient  form  for  computation  by  means 
of  desk  calculators,  was  presented  by  Mr.  Emanuel  M.  Sodano  at  the 
Xlth  General  Assembly  of  the  International  Association  of  Geodesy  and 
Geophysics  in  Toronto,  Canada  in  1957-  The  results  represent  the 
gradual  extension  and  accumulated  improvements  of  the  original  Army 
Map  Service  Technical  Report  No.  7. 

This  modification  contains  a  more  stable  formula  for  azimuths  and 
an  alternative  formula  for  very  short  lines.  More  general  and  accu¬ 
rate  formulae  for  both  long  and  short  lines  are  given  herein  than  are 
contained  in  Technical  Report  No.  7*  The  complete  theoretical  deriva¬ 
tion  starting  with  a  rigorous  modification  of  Helmert's(l)  classical 
formulas  are  given.  The  final  non-iterative  formulas  have  been  ex¬ 
tended  through  terms  equivalent  to  the  second,  fourth  and  sixth  powers 
of  the  eccentricity  of  the  spheroid,  and  therefore,  may  be  shortened 
according  to  the  required  accuracy. 
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The  solution,  which  requires  no  special  purpose  tables,  is 
accurate  to  at  least  the  tenth  de'*-’»»l  place  of  radians  for  +''r' 
azimuths  and  the  arc  distance*  If  the  final  formulas  are  shortened 
to  the  second  and  fourth  powers  of  the  eccentricity  respectively, 
the  results  are  accurate  to  seven  and  nine  decimal  places  of  radi¬ 
ans  respectively,  even  for  distances  circumscribing  the  earth. 

In  Section  IV  the  formulas  for  the  solution  of  the  Inverse 
Geodetic  Problem  have  been  adapted  to  electronic  computers.  These 
formulas  were  derived  from  the  basic  formulas  of  Section  III.  A 
solution  of  the  Direct  Geodetic  problem  is  given  in  Section  V. 

The  formulas  are  adapted  to  electronic  computers. 


SECTION  ni. 


A  RIGOROUS  NON- ITERATIVE  PROCEDURE  FOR  RAPID  INVERSE  SOLUTION 

<3> 

OF  VERY  LONG  GEODESICS 

2.  Preliminary  Modification  of  Helmert's  Iterative  Solution 


e 

e» 


fix 


and 


eccentricity  of  the  spheroid  * 
second  eccentricity 


V 


ao  -  ho 


n — 77 

a0  -  bo 


2*  a 


"bfT 


semi-minor  axis 

absolute  difference  of  longitude  on  the  spheroid, 
between  the  given  endnoints  of  the  geodesic, 
parametric  (or  reduced)  latitude  of  the  westward 


'•ni  eastward  ond-oint.- ,  respectively. 
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The  relationship  between  parametric  latitude  and  geodetic 

latitude  is  given  by  the  equation  tan^  -  Us  B(l-f)  '-here  t 

is  the  spheroidal  flattening. 

A  -  difference  of  longitude  (approximately  L)  on  the 

reduced  sphere,  for  which  a  progressively  better 

value  is  found  with  each  repetition  of  the  follow- 


cos  /0 
sin  & 

fo 

sin  2§0 
sin 

OOM/SQ 

sin  2fi0 

cos  20- 
cos  ha- 
cos  6  a- 
A  * 


inp  iteration  process* 

sin/3^  sin  A  ♦  C0SA  tfosA  cos  X 

(sign  of  sin  A)  V 1  -  cos2  §0 


positive  radians 

2  sin  /0  cos  j0 

3  sin  /0  -  U  sin3  j0 

(cos  A  cos  /&2  sin  X)  *  sin 
1  -  cos2/60 

(2  sinA  sin  A  7  si*1*  A)  -  008  /o 
-1  +  2  cos2  2  c- 


U  cos3  2<r  -  3  cos  2  ©" 


-  42cin^0  *  ■‘"V. 


B--  bin2/?  -  2?£^  .in1* /?„ 

16  ^ 0  32 

fi2  |J|  h  yn 

C*  ■  g„g  sin  fio 

1ZT  r 

T  -  A»/0  -  B'  sin /0  cos  2«r  ♦  C'  sin  2jf0  cos  U  tr 
Next  approximation  to  A  •  [(L  ♦  T  cos  A)  radians. J 

After  a  sufficiently  accurate  X  is  found,  and  using  the  set 
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of  valuer  from  the  last  iteration,  the  geodetic  distance  (S)  and 
azimuths  (or)  between  the  endnointrare  obtained  as  follows: 

Ao  =  1  +  ^2-  sin2^0  -  3e^__  +  |e*  sin6^o 


—  sin2/J  -  e'U  sinM  ♦  tfe’6  sin6^ 

It  /  0  TT  0  “H? 

S5sinU/^°-  ^  sin6^° 

e'6  _,_6>3 


no  * _ .  sin  A 

1‘-  36 


f  S  -  b0(A0/0  ♦  P0  sin  /0  cos  2cr  -  C0  sin  2jf0  cos  b  er 

♦  D0  sin  3$)  coa  6®“) 

1  f  ^1  2  =  tan  ^2  cos  /^-cos  A  sin  ' 


sin  * 


k/ }  a  ni  ?  COS.A-  cos  A  ban/^ 

sin  A  -I 

£*]_;>  ir-d  ranqe  from  0°  to  180°  and  180°  to  360°, 


:*•  ;  eotj vely,  clockwise  from  north. 


L^t  it  be  assumed  that  the  true  value  ofAis  known  (that  is, 

V value  that  would  result  from  an  i  nfinite  number  of  Helmert  ap¬ 
prox!  nations)  and  let  this  true  value  bo  reoresented  by  the  piven 
.•ins -lute  difference  of  longitude  on  the  spheroid  plus  a  quantity  x 
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which  will  be  determined  later. 

Thus:  A  -  (L  +  x)  . 

It  will  be  evident,  later,  that  x  is  a  very  small  positive  quantity 
of  the  order  of  e^,  and  therefore  veil  suited  for  setting  up  a 
convergent  nower  series  in  x  for  each  expression  contained  in  the 
Helmert  procedure.  For  example,  from  the  above  assumed  equation, 
the  following  is  derived: 

cos  A  ■  cos  (L  ♦  x) 

■  cos  L  cos  x  -  sin  L  sin  x 


-  (sin  L)  (x  -  ...) 

"•herefore:  cos  A  ■  (cos  L)  -  (sin  L)  x  -  £(cos  L)  x2  +  .  .  . 

There  is  thus  available,  at  the  outset,  a  series  for  the  true 
cos  A  with  which  to  begin  the  Helmert  solution  and  develop  it  in 
rover  series  In  x  in  its  entirety.  The  process  consists  of  sub¬ 
stituting  each  new  series  into  the  succeeding  Helmert  expressions 
as  reouired.  For  convenience,  the  following  additional  notation 
will  be  used: 

N  •  e '  £  (e '  *  e) 
a  ■*  sin^i  sin/^j 
b  *  cos/3}  cos  $2 
cos  $  *  a  ♦  b  cos  L 
c  *  b  sin  L  esc  § 
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m  =  1  -  e2 

h  a  e,2m 

F  *  n  cct  §  -  a  esc  § 

T \  -  (tan /5?  cos/^i  -  cos  L  sin/^)  *  sin  L 

"■?  *  (sin/^2  cos  L  -  cos/^2  tan^i)  |  sin  L 

listed,  bolov*,  in  the  same  sequence  as  the  corresponding  Helmert 
expressions,  is  the  complete  set  of  series  through  Tcos/^#*  The 
extent  of  the  movers  of  x  is  such  as  to  permit  accuracies  of  the 
e*  order  in  A  ,  for  subsequent  application  to  the  distance  and  azi¬ 
muths  to  the  sane  decree  of  accuracy  as  the  reference  Helmert  iter- 


.*  t::  or.  form. 


on: 


/, 


>  r 

h 


- 

sir;  A 

■  p,  „ 


C'  -;  ?  <r 


-  (cos  fi)  -  (c  sir.  $)  x  -  l(c-  cos  jf  +  P  sin  J?)  x2 

-  (rin  /)  *  (c  cos  jjT)x  -  f(c~  sin  jf  -  P  cos  /)  x2 
*  /  +  (f0-  §')*§+  arc  sin  ["sin  (^0  - 

=  jf  ♦  arc  sin  ["sir.  cos  -  cos  $5  sin  J  J 
=  f  *  (c)  X  *  s(P)  X2 
=  2  sin  <f0 cos /<, 


=  (sin  L)  *■  (cos  L)  x  -  sin  L)  x2 
=  (o)  +  (°)  x  -  f(cm  ♦  3=  P  cot  $)  x2 
=  m  -  (2cF;)  x 

r  A  (m  cos  -  2P  sin  (cm2  sin  0"  +  h  ctnF  cos  jT 

n  m1-  __  “ 

-  ucP^  sin  jf)  x 
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cos  I4  or 

1 

«c  _ 

-r2 

(m2  -  2m2  sin2  £  -  8mP  sin  JT  cos 

+  8P2sin2  JT) 

A* 

e2 

(12®  -  8h 

+  3h2)  +  (.e*2  cP)  x 

128 

8 

B' 

c  _ 

(2h  -  h2)  - 

®i(e'2  cP)  x 

32 

T 

C' 

e? 

ac 

_  (h2) 

256 

A'/o 

2  2 
«  l—  (12®  /  -  8h  0~  +  3h2  /)  ♦  f_  (lfiKc  -  he 
.128  16 

+  2iicY§) 

+ 

(NP)  X2 

2 

^  e2 

-B'  sin  Ceos  2  ®"  =  — 

(-  8h  sin  §  cos  $  ♦ 

I6e  '2  P  sin2  § 

—  w 

128 

*  Uh2  sin  £  cos  $  - 

86'2  hP 

sin2  £) 

-  (he)  x 

16 

C '  sin  2jT0 

e2 

coslr-  -8 

(h2  sin  £  cos  £  -  2h 

2  sin-^  ^ 

cos  § 

-  8e'2  hP  sin2  jT cos2  ^ 

+  Se'^P2 

sin^  £  cc 

Tco4 

=  J-_  (l28Nc/- 

1  ?P.  ^ 

8hc0"  -  8hc  sin  JT  cos 

/  «■  l6e  ’ 

2cP  sin2 

3h2c^T  ♦  5h2c 

sin  §  cos  jT  -  2h2c  sin^  §  cos 

■•i 

- 

8e ,2hcP  sinf 

_0"-8d2! icP  sin2  0  cos2  0 

+ 

Pe'^cP2  sir.^ 

0" cos  ^  (}6Nc2 

+  l6NPjf 

-  2hc2 

- 

hP/  +  2e,2c2 

PjT  -  hP  sin  £  cos  J0" 

«•  2e ,2P2 

sin2  > 

- 

2-  (Ncnjf  -  3PcP  +  3’!cD/  cot  0)x2 

L •  Derivation  cf  the  Unknown  Quantity  X 

Since  the  substitution  into  the  Helmert  iteration  began  wii!  an 
alrebmic  series  representing  the  true  X,  the  next  approximation  tc 
A  must  of  necessity  be  its  equalj  that  is: 

The  next  approximation  to  A  *  the  starting  true  ^ 
or  L  +  Y  cos  SQ  B  L  +  x 

and  therefore  *£  cosS0  *  x. 

By  replacing  T  cos /3Q  with  its  corresponding  power  series,  the  above 
equation  takes  the  following  quadratic  form: 

Ql  +  Q2x  *  Q3X?  “  x 

for  which  the  reouired  solution  o  f  x  to  the  proper  order  is 

x  *  Q;j_(l  +  Q2  +  Q|  ♦  QiQ 3). 

Finally,  substituting  for  Qi,  Q2  and  Q3,  produces  the  following  end 
result : 

x  =  1_£  |^12  Rh/  +  l?t'e?-r.2c?f~  Bh/  -  Ph  sun/  cos/  ♦  12Be?K2F/’ 

+  16e,2P  sin2/  +  ISfieMc/ ~  2he2Nhc2/  +  3h2/ 

-  8e2Nhc2  sin  $  cos  /  +  £h2  sin  /  cos  /  -  2^5 

-  2tv  sin-/  cos  £  +  (16e2e,2P  +  liLPeWi^)  c2P  / 

-  16e2PhP  /  +  l6e2e,2Nc2P  sin2  /  -  8e,2hP  sir.2  £ 

16e2NhP  /  sin  /  cos  /  -  19Pe^N^c2P  jjP  cot  / 

-  8e,2hPsin2/ cos2  /  +  128eV^p/^ 

3  ?e2e  ,2PP2  £  sin2/  +  Be  '^p2  sin-^  /  cos  jfl 


IP 


The  above  rigorously  developed  e xpression  is  completely  non-iter¬ 
ative,  since  it  requires  only  the  Riven  snhero:idal  longitude.  It 
therefore  Permits  a  direct  evaluation  of  the  ultimately  true  A  (that 
is,  L  +  x),  extended  in  this  case  through  terms  equivalent  to  the 
e^,  and  e^1  order  consecutively,  in  accordance  to  the  accuracy 
that  may  be  desired.  Furt' ermore,  it  represents  the  algebraic  solu¬ 
tion  of  the  hitherto  unkovm  quantity  x  used  in  the  power  series 
version  of  each  of  the  intermediate  Helmert  expressions* 

5.  Determination  of  Geodetic  Distance  and  Azimuths 

The  non-iterative  expression  that  has  been  developed  for  x 
suggests  at  once  a  numerical  solution  of  distance  rnd  azimuths  where¬ 
in,  using  the  resulting  true  value  of  A,  only  a  single  evaluation  of 
Helmert ’s  original  formulas  is  necessary.  An  illustrative  example 
by  such  a  procedure  is  given  in  po 

On  the  other  hand,  instead  of  reverting  to  functions  of  the 
true  A,  the  distance  and  azimuths  themselves  can  be  expanded  non- 
iteral  ively  into  rover  series  of  x  with  coefficients  in  terms  of 
the  riven  spheroidal  difference  of  longitude.  This  is  accomplished 
below,  but.  limited  to  the  order  of  accuracy,  since  this  manner 

of  obtainir.r  the  distance  and  azimuths  through  would  require 

each  component.  series  to  one  hirher  rover  of  x  than  was  necessary 
f or  A  .  Analn,  the  series  are  developed  in  the  same  sequence  as 
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tho  corresponding  Helmert.  expression. 


A0  r  ^  (6U  ♦  l6h  -  3h2)  -  i  (e,2cP)  x 

Ec  =  L  (hh  -  h2)  -  \  (e,2cP)  x 
1c 

h2 

C°  ’  128 

A0  f  -  i-  (61*/  ♦  I6h/ -  3h2/  +  w  (Lc  ♦  he  -  2e«2cP/}  x 
6i* 

♦  i  (p)  x2 

Pcrin/,ccs  2<r»  ~  (l6h  sin  / cos  /  -  3?e*2P  sin2/-  lah2sin  /  cos  / 
+  8e ,2hP  sin2  /  ♦  \  (he)  x 


-C(  rir.p/crr  1^=  ^  (-h?sir/  cos /  ♦  ?  r''Viri P«- ,2V^  ?:in2/  c-s2/ 

-Pe  *Lp2  sir./  cos  /) 

'  b  ' 

r  »-  (61/  ♦  l6h/  ♦  l6h  sin  /  cos  /  -  22n  '2P  pin2/ 

oli  * 

-3h?/  -  5h?oln /Vv/  •  'h2:- in/cr-r  / 

+  P<»*2hP  sin2/  +  8e ‘2hl;  sitr/ccs2/ 

-  ee,ijP2  sin/  cos/)  *  (2c  +  he  -  e  ,?cP/)  x 

♦  b°  (P)  x2 

T  / 


Up  cos/^ 


T1  1  sin  L  cos^J  p 


\  /  T'i  Up  cos  L  cos/^  \  ;1 

?. J  ^2  sin2  L  2  sir,2  1,  cos /&J 


cot«Xp_^  ■  Up 


/Ul  eos^  \  J  U 
(  sin  L  cos^  I2  si: 


COS  L  CCS // 


in2  L 


2  j  X2 


2  sin2  L  cos,^ 

The  x  and  xc  for  the  above  formulas  of  distance  and  azimuths 
can  be  substituted  either  numerically  or  algebraically  using,  in  this 
case,  only  the  first  6  terms  of  x  for  accuracies  equivalent  to  the 
order.  The  algebraic  substitution  gives  the  following  final  ex¬ 
pressions  : 

S  *  W  *  61je2Ne2/  +  l6h/"  ♦  l6h  sin  /  cos  / 

-  32e,2P  sin2/  +  6Ue^,2ci/'“  3h2/  ♦  (32e2N  -  Ue2)hc 2jf 

-  lie2hc2  sin  /"  cos  /  -  S>h2  sin  /"  cos  /  ♦  2h2  sin3jf  cos/" 

+  (96e^N2  -  3?o?e  ,2N)c2P  /^+  8e2e,2o2P  sin2 ^/ 

+  f-e,2hP  sin2  +  6e,2hP  sin2 /"  cos'/ 

-  fle'kp2  sin^  /  cos  /J 


cot^. 


Ul 


e2Nc/?2  cos/^ 
sin  L  cos,$ 


! ,  o  i 

e“N‘-e-'/l!2  cos/-?i 


2  sin  L  cos^2 

e2hc/j?  co^[  +  e2hcU2  sin/  cos/"  cos/(^ 

16  sin  L  cos/£  16  sin  L  cos^ 

t  N2cl/®!l2COE^  «,  e2e,2cFU2  sin2/  cos/j^ 
sin  L  cos  /?2  S  sin  L  cos 

e^N"c2/^L'2  cos  L  cos/^  e^2c2/H'2 


2  sin2  L  cos,^  2  sin2  1 

The  corresoondinr  cotoi'2-l  Is  obtainable  from  the  above  by  inter¬ 
changing  tT^  with  Up  and  with 

IS 


Thu?,  progress i vely ,  there  have  been  (level ored  three  rigorous 
methods  for  determining  geodetic  distance  and  azimuths  non-iters  ti vely: 
as  a  function  of  the  true  A,  as  a  newer  series  in  x,  and  culminated 
by  an  explicit  expression  in  essentially  the  given  snheroidal  lati¬ 
tude  and  longitude  of  the  endpoints.  For  shorter  lines,  or  for  re¬ 
duced  accuracy  on  long  lines,  terms  may  be  still  further  eliminated 
according  to  the  next  higher  covers  of  e2,  e’2,  'n  and  x,  or  equiva¬ 
lent  combinations  thereof. 

6.  Other  Won-Iterative  Solutions 

The  distance  and  azimuths  by  the  original  Helmyrt  method  are 
essentially  functions  of  elements  in  the  following  spherical  triangle: 


/. 


where  ^ 

§o  "  T2  ~  T1 
2  <J~  ■=  T 2  +  ,7\ 

and  .Ap,  Aj,  and  0"->  are  regarded  as  negative  or  positive  accord- 
i r.«*  to  whether  they  are  west  or  east  of  the  perpendicular  arc  90r -/3L. 
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(^or  this  specific  configuration,  therefore,  A  arri^  actually  rep¬ 
resent  the  sum  of  the  absolute  con^onents  and  2  o'  t  he  di  fference.) 

Since  Helmert's  method  of  successive  approximations  can  only 
determine  A  first,  the  subsequent  solution  of  the  above  spherical 
triangle  would  always  begin  with  A  and  the  known  and  ^ .  The 
present  paper,  however,  has  developed  net  only  a  non-iterative  ex¬ 
pression  for  A,  but  also  independent  never  series  for  the  various 
elements  of  this  spherical  triangle  or  functions  thereof.  Therefore 
the  combination  of  ways  to  compute  quantities  leading  to  the  dis¬ 
tance  and  azimuths  is  increased  considerably.  In  addition,  the  x 
for  such  series  can  be  substituted  either  numerically  or  algebraic¬ 
ally,  in  the  manner  shown  fer  the  distance  and  azimuth  series  in 


Tie  above  potentiality  for  increasing  the  number  cf  non-itera¬ 
tive  solutions  may  be  seen  from  the  expressions  (1)  belcv:,  where-ir* 
the  x  and  7?  of  the  j?o  series  are  algebraically  eliminated. 

|/0  .  /  ♦  (C2/)  ♦ 


16 

e?e'2  f 


16 


c2sin0cosjO  -  e  e  cksinj^cos^) 
16  16 


(1) 


e^e  /r\  + 

"  (airsinf)  +  - 


2  (e?df)  -  (c^cot/)} 

Thie  computed  value  of  is  then  combined  with  and  ^ 
to  obtain  's,  followed  by  /S0)  ?0F',  A0,  B0,  C0  and  finally  the 
geodetic  distance.  When  adopting  such  varied  procedures  for  solving 
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the  reference  triangle,  care  should  bo  taken  to  avoid  ton .ulati one 
vhich  lead  to  a  veak  determination  of  required  quantities.  There 
difficulties  may  most  likely  occur  at  extremes  of  latitude,  loi  ,i~ 
tude,  or  azimuth. 

The  non-iterative  series,  too,  are  functions  of  element:  of 
a  spherical  triangle,  but  defined  by  and  {3?  and  the  gi  ven 
longitude  L.  This  amounts  simply  to  a  substitution  of  L  for  A, 
which  results  in  a  spherical  triangle  with  parts  corresponding  as 
follows: 

Series  ^1  @2  L  /  c  Up  U2 

i’.elmert  /ip  /?2  A  f0  cos  /?0  cot°(p_p  cpV*p_p 

Starting  wj  th  the  given  ana  L,  the  values  of  all  quanti¬ 

ties  used  in  the  non-iterative  series  may  thus  be  solved  trigono¬ 
metrically  in  various  orders. 

Tt  is  also  to  be  noted  that  in  the  relrtior.  A  =  (L  +  x),  if 
y  is  assumed  to  be  zero,  L  is  cnnsider°d  to  be  equal  to  A  . 
Therefore  in  the  various  mover  scries  in  x,  the  constant,  term 
can  ror>res*'nt  the  true  value  cf  the  series  by  simply  replacing 
functions  of  L  with  A.  This  is  well  illustrated  by  the  A0 
series  in  x  in  -•ar'Trnnb  :  and  itc  counterpart  in  :v*ra  \;v.:  ~ 

■  •  r  nv'er.  The  principle  can  veil  be  incorporated 

un  com} -utat j or,  fonts,  such  as  tuo  one  cn  the  next  ■_ \: ■.  ]  *■  ' r  ' . 
apolied  to  H,  I,  etc. 
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7.  Numerical  Illustration  of  a  Sample  Solution  (International  Ppherold) 


'L  •>  absolute  difference  of  lonritude 
Given  ^  Pi  c  latitude  of  westward  point 
>  latitude  of  eastward  point 
=  0.99663  29966  tan  Pi 
tan (2 2  ~  0,99663  29 96'6  tan  B2 
js/?i  •  1  •  t  ■\/l  +  tan'v^ 


cos 

co; 

sin 


/?2  1  -  ♦  V1  ♦  tan^ 

/^1  =  tan  £%\  cos^i 

sin  $2  =  tan  ^2 

a  *  sin>5}  sin 

t  =  cos/^i  cos 


sin  L 


106° 
20°  N 
145°  N 

0.36275  l*7L53 

0.99663  29966 
0.95006  23275 
0.70629  61969 
0.35100  26695 
0.7059]  33515 
0.2u071  63353 
0.665 65  55515 
0.96126  16959 


cos  I 


fo. 27563  73558 


cos  jf 
rin  jf 


1 


A 


C 


a  +  b  cos  L 

(sipr;  of  sin  L)  -  cos^jf 
positive  radians 
(b  sin  L)  t  sin 
A? 


0.05716  6736i3 
0.99836  36996 
1.5135?  63766 
0.61,109  99269 
C. y 1100  91163 


{cos  j7"  -  (cos  pj  -  206.59 


0.00675  li?028 


D  «=  -a  (O..U0106  12630) 
=  -a (6. 79950  93666) 
F  -  (3.9665  20659)C 
O  c  ^  ^  sin 


0.09656  7615? 
■0.19251  65307 
0.02692  7762? 
2.2956?  57386 
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'TPO 


I A  [/( 2  37. 2388918  ♦  B)+sin/(C+D) 
♦0(F+E)]j  7  70519. 51115 


0.00326  58 167 


A  -  L  ♦  x 

“in  A 

cos  A 

cos  $r  °  a  ♦  b  cos  A 
sin  </r  ~  (sign  of  sin  A)  v/T  -  cos  % 

=  positive  radians 
in  -'t.  -  (sin  /Q  cos/0)  i  0.5 

eos/y  ~  (b  sin  A)  *  sin 

'1  *  1  -  COS2/^0 

;;<r  -  (2a  -  q  cos  $0  *  q 
i.i  <r=  (cos2  2<r  -  0.5)  f  0.5 

H  **  6356911. 9W  +  10756 .165q  - 

■3/50o? 
i  «  10756. I65q  -  lfi.200q2 
•T  =  2.?75q2 

d-  j_5  =  (tanjjj  cos /3\  -  cosAsin^i) 

sUTA 

=  (sin cos  A  -  cos/?2t-an/?i) 

sSrTTf 

=  Clockvise  from  North,  in  quad 
I  or  II  for  cot  ♦  or  -,  reap. 


(A  - ' 


a 


1C6°11 '13'' .62305 
0.96035  63907 
-0.27677  510U6 
0.05509  ?U?83 
0.9 96 UP  09630 
1.51567  00535 
0 .11002  ?ii687 
0.6U0U?  0753? 
0.58986  12195 
0.76108  80.731 
0.15651  2697? 

6  36  32 5 1.8  hi 
6338.31? 

0.79? 

1.07U55  96397 
-O.U72U5  22891 
U2°56’30,,.03557 


or 


rv'1  ■:  V: 


=  Clockwise  from  North,  in  quad 

ITT  rr  IV  for  cot  ♦  or  -,  rear.  295°17  '18” .59865 


■  K/o  *  I  sin^j  cos2  - 


96U9U12 ,85Lm. 


s;  r!  ?jj£.  C03h  <T 
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8.  Numerical  Coefficients  For  Other  Spheroids 

The  illustrative  solution  piven  in  the  preceding  section 
contains  fixed  numerical  coefficients  which  are  functions  solely 
of  the  size  and  shape  of  the  International  spheroid.  The  alpebraic 
exoressions  of  these  coefficients,  together  with  their  values,  are 
shmn  below  in  the  order  of  appearance  in  the  samle  solution. 

For  any  other  spheroid,  the so  expressions  can  be  quickly  re-evalu¬ 
ated  once  and  for  all  and  substituted  for  the  corresponding 
International  values.  (Note:  e2N  s  flattening.) 


♦  Vl  -  e2 

■ 

0.99663 

29966 

(l6e2N2  +  e'2)  i  e‘2 

- 

h.9866 

206U9 

2e»2  i  (l6e2N2  +  e'2) 

s 

0.10108 

12630 

l6e2N2  a  (I6e2h2  ♦  e’2) 

■ 

0.799L5 

93686 

16o2N2  i  e’2 

a 

3.9865 

206h9 

(16N  -  e'2)  i  (I6e2!:2  +  e'2) 

a 

237.2386 

918 

16  i  e2(l6e2N2  +  e'2) 

a 

7C519.511U5 

*>o 

a 

6356911.91:6 

bce'2  *  li 

= 

10756.165 

3b0e'l  i  61 

r 

13.650 

b0e  i  16 

a 

18.200 

b0e’k  1  128 

a 

2.275 

9 •  Additional  Notes  on  Computational  Procedures 

Although  the  illustrative  solution  piven  in  ;  — . .  ;  7 
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is  primarily  intended  for  accuracy  equivalent  to  the  order,  it 
easily  lends  itself  to  any  required  degree.  This  is  accomplished 
simply  by  adding  or  subtracting  appropriate  terrna  of  x,  H,  1,  J, 
and  S.  The  extended  terms  are  given  in  the  latter  part  of  Section 
III,  ".-rorr"''}:?  );  and  2  respectively.  For  short  lines  or  re¬ 
duced  accuracy  on  long  lines,  x  on  the  International  spheroid  be¬ 
comes  merely  (kjf  $  297)  and  all  terms  in  q^  are  omitted,  with  the 
consequent  elimination  of  many  other  supporting  quantities.  Simi¬ 
lar  savings  are  realized  for  other  forms  of  solutions  presented 
herein . 

For  short  lines,  the  resulting  small  §  is  computed  more 
accurately  from  sinjT  obtained  as  follows: 


Similarly,  sin  §Q  is  obtained  as  above  by  replacing  §  with  £0  and 
L  with  )\ .  .si  either  case,  squaring  the  small  sines  under  the 
radical  increases  their  significant  decimal  places. 

If  the  numerator  of  x  is  to  be  cumulated  in  a  ter.  digit 
calculator,  9  decimal  places  should  be  allotted  to  §,  sin  §  snd  0, 
but  only  7  decimals  to  their  multipliers.  However,  when  the  value 
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of  0  is  10  or  greater,  decrease  its  decimal  places  accordingly  and 
increase  those  of  F  and  E  correspondingly.  For  a  smaller  calcula¬ 
tor,  reduce  all  decimals  equally. 

Use  co-function  of  tan/9  or  cot  of  when  their  values  are  too 

large. 

lhus  cot  /^n  -  — cot  Bn  tan  c<  ■  — - - 

+  -/  cot* 

The  accuracy  of  geodetic  distances  computed  through  the 
e^,  and  e^  order  for  very  long  geodesics  is  within  a  few  meters, 
centimeters  and  tenths  of  millimeters  respectively.  Azimuths  are 
good  to  tenths,  thousandths,  and  hundred  thousandths  of  a  second. 
Further  improvement  of  results  occurs  for  shorter  lines. 

■  Some  of  the  terms  in  the  sample  solution  of  paragraph  7 
have  been  grouped  for  ease  of  computing  by  desk  calculator.  For 
electronic  computers,  however,  the  terms  are  best  left  in  series 
form,  thus  being  ideally  suited  to  adding  or  removing  them  accord¬ 
ing  to  accuracy  requirements. 

10.  Antipodal  Points 

In  the  various  series  that  have  been  presented,  $  represents  a 
spherical  arc  distance  which  varies  from  0°  to  180°  and  even  to 
360°  according  to  whether  the  geodetic  line  is  very  short,  half 
around  the  earth  or  completely  around  it.  At  these  specific  instan¬ 
ces,  quantities  such  as  esc  jf,  cot  and  P  approach  infinity.  For 
the  case  of  the  very  short  lines,  this  condition  is  equalized 
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* 

where,  as  noted,  the  substitution  of  the  ^cos^0  series  given  at 
the  end  of  paragraph  3  will  now  be  in  terms  of  I*,  and  <  instead 
of  L  and  x.  Solving  the  above  equation  for  Z  (this  time  through 
only  the  order  of  accuracy)  gives: 

16(1^1^,)  +  (16e2Nc/  -  e2hc/  -  e^hc  sin/  cos/  +  2e^e'^cP  sin^/)n 
16(1  -  e2Nc2  -  o2Nl/)n 

where  the  subscripts  n  to  the  parenthesis  indicate  that  c,  /h,  P, 
etc.  are  functions  of  Ln  instead  of  L«  This  time,  the  denominator 
of  the  expression  cannot  be  algebraically  divided  into  the  numerator, 
because  the  e^NP/  term  is  relatively  large  for  nearly  antipodal  lines 
With  the  above  correction  x  to  an  arbitrary  but  sufficiently 
accurate  value  Lp,  the  true  A  of  antipodal  lines  is  essentially  ob¬ 
tained  again  non-iteratively,  and  therefore  more  rapidly  than  by 
numerous  individual  successive  approximations.  Thus,  also,  a  pre¬ 
vious  it"  longitude  discrepancy  noted  by  Mr.  H.  F.  Rainsfon/^for  a 
line  of  about  l?9°Li6 '18"  longitude  would  be  resolved,  3h  this  con¬ 
nection,  appreciation  is  expressed  to  Mr.  Rainsford  for  his  interest 
in  the  subject  which  resulted  in  profitable  corresnondence. 

SECTION  IV 

TAPULAR  AND  ELECTRONIC  COMPUTER  ME  TROD  FOR  NON-ITERATIVE 
SOLUTION  OF  GEODETIC  INVERSE,  RASED  ON  LGMNO 'S  PAPER 

Due  to  their  series-like  nature,  the  formulas  given  in  this  sect-5 
for  distance  and  azimuth  are  more  adaptable  to  electronic  comuier 
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programming  then  the  corresponding  closed formulas  of  Section  in. 

Thir  r.etl.ed  (unlike  the  one  just  discussed)  doer  not  l.-'.vu  1 
restriction  that  P]_  and  L]_  must  be  the  latitude  arid  longitud.  ,  r._- 
pectively,  of  the  westward  point.  Here,  R]_  and  are  the  gee  r:n-\hic 
latitude  end  longitude,  respect! vely,  of  any  point. 

The  distance  equation  of  this  section  was.  derived  by  making  the 
* o2 3 ovj np  fibrti tutions  into  the  distance  S  equation  on  pope  Ui. 

f  =  1 

e'-n  =  h  ('-'here  e1^  was  expressed  in  terras  of  f)  ?  (1) 

m  =  1  -  C"  J 

P  =  (1  -  c^)  cot  §  -  a  esc  £ 

The  expression  ( A  -  L)  of  this  method  is  equivalent  to  "x"  on 

rr-  •*c-  I  •  The  series  for  (X  -  1)  was  derived  by  nnkin  -  the  suhsti  tu- 
tions  (!)  into  the  equation  for  "x1'  on  p?  e  11.  The  corf uinta on 
t’.-r  :  fer  -this  .-.othed  is  as  follows : 

?] ,  L]  -  reoyrarhic  latitude  and  Ion --i  turn,  rerocct:'  '■  ly, 
of  any  noint- 

Bp,  Ip  =  geographic  latitude  and  longitude,  respectively, 

of  ary  other  noint,  inn-antipodal 
Lati  t.r’er  and  longitude?  considered  ( - )  a  cth  an  1  east,  (-)  south 
&rvi  vTC5s^ 

I.eo.iiiu.-ds  Of  ,  3  =  azixutb.s  clockvi  re  from  north  and  distance 

between  roints,  respect! vely. 
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wherej  ap,  b0  *  semi-major  and  semi-minor  axes,  respectively,  of 
spheroid 

f  ■  spheroidal  flattening  ■  ^1 
P  *  number  of  seconds  in  one  radian  *  206,?£b»C0("'2i> 

L  *>  (Lj-Li)  or  (Lg-Li)  *  [sign  opposite  of  (L?-Li)  J  360° 

"i  i  ■••'rr  L  has  an  absolute  value  <  or  >160° j  according  to 
whether  the  shorter  or  longer  geodetic  arc  is  required,  However, 
for  meridional  arcs  (|l|  *  0°  or  180°  or  360°)  use  either  L  but 
consider*  it  as  (  +  )  for  the  shorter  and  (-)  for  the  longer  arc. 

tun jS  M  tan  B  (1-f)  vhenjB|£  L£°  or  cot^  ■  when 

|r^>  \6° 

o  *>  sin  fi  i  sin 
V  «  cos  /}\  cos/?2 
cos  jT  -  a  ♦  b  cos  L 

sir.  /  -  *  >/(sin  L  c°^)2  +  (sin^  cos^  -  sin^  cos/^>  cos  L)*- 


The  sign  o'"  sin  jp  is  (*)  or  (-)  according  to  whether  the  shorter 
cr  the  longer  arc  is  required.  'The  quantity  under  the  radical  and 
its  root  must  be  onmou-ed  by  floating  decimal  to  obtain  sin  §  to 
full  accuracy  for  short  lines. 

•f  -  positive  radians  (obtain  reference  angle  from  sin  § 
or  cos  §  .  '  has  smaller  absolute  value.  ) 

c  =  (b  sir  L)  7  sir.  § 
r.  =  "  -  .■ 


2d 


cot*<]L_2  =  (tan/? 2  -  cc 

cot  <^2-1  “  (cos  ^£00^2  ~  tan/ 

If  |  cot  o<  |  >  1,  divide  result  into 


mm 

L 

Sipn  of 
tan  or  (cot) 
<X  1-? 

+ 

+ 

V 

•i. 

T  - 
j.  J. 

- 

+ 

Ill 

- 

IV 

7'>r 

’■nr  *  ~  » .  r-  ].  C.TCF-  ('D*.*"1!’  the 

r.ujier-1  ten 

of  ent„  C< , 

•I  roferer  c« 

1 1 .  Ex+  ar. fii on  and  * 

*r>-  r 

Mel 

•c-*t  for  don-It; 

.v;.  ••■.  :  ■  lor.: 

I'T.re?  fieri  Decl 

'•cccrtov  ir.  i 

»r.  _ 

t  ■  <  V.fi  1  :>cc  r .. 

..-  Yu;  • 

o  ■■■'■?  -  'i 

!  g  1-  Y 

t-  do  :i  '■ 

F- 

T  '  finrie?  liitr.  ■.•nc--.. 

:• : 
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+  m 


T 


3  an  (  1  IT  \ 

F0 

+  a  j^(f+f2tf3)sin/  ♦  (-’f?-f3)/2cse/  *  ff3/3csc/ cot/j 

f  Bf-if?-if3)/  +  (-^f-^f?-?1f3)Pin^r  cos/ 

+  (sf2*’f ^J^cot/  -  3 /3  _  i-f/*Cot2/j 

+  a‘_  |"(-?f2-f3)sin/  cos/"  +  |f3/3  esc2/  ♦  £f3/l 

+  m2f  (+  ^f2  •*•  3  f3)/  +  (  3^f2  +  3f3)sin/  cos/ 

+  (~|f2  -  ^  f3)/cot/  +  (-  ^  f2  -  if3) sin/ cos3/ 

♦  if3/  cos  2/  ♦  .^f3/3  +  ^  f3/  3cot2/j 

(|f2  ♦  ^f3)/2  esc/  +  (|f2  +  f3)sin/  cos2/ 

-  3  f3/  cos/  -  2f3/  esc/  cot/| 

-in  -  if3/-  |f3sin/  cos/  -  if^Poac2/  +  f3sin3/  cos/| 

O*  -  J 

-  3  f3/2  esc/  +  if3sin/  cos  2/+  3f3/  cos/ 

U  4 

>■  « 

+  3  f3/3  esc/  cot/  -  ^f3sin/  +  if3  sin"/ 

2  ‘ 

^  ^  f3/  «•  3  f 3/2  cot/  -  ^  f  3  sin/  cos/ 


+  am 


+  a* 


+  am 


+  ^  f3sin/  cos3/  -  ^f3/  cos2/  -  3  f3/1  -  f3/3cot2j' 


[* 


+  A  f3sin3/  cos3/ 


1? 


+  a3  |  -|-f3  sin/  -  2  f3  sin3/ 

a 
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The  maximum  values  for  the  f3  term  of  have  been  found  for 
various  lengths  of  arc  /  and  are  recorded  Id  Appendix  I,  Part  A. 

3ii:ilsrl.v,  the  scourac.v  of  the  ^  series  may  be  extended  by 
adiinr  the  f-*  term.  The  series  t1  en  becomes : 

*  ^(f+r+f3)  /j  *  a  |'(~;1:f?-f3)sip/  +  (-f?-lf3)/?  esc/ 

+  2  f3.Pc*cf  cot/ j  +  m  £(-  |  f?-3f3)/  +  (£f2+i,f3)sir/cos/ 

♦  cot£  -  ff3/3  -  3  f  3/3Cct//j 

+  ra?  ^  TP  fV  “  x/3  sin/  cos/  +  if/3  -  ^f3si.n3/  cos/ 

-  ffi/ent/  +  ;'f3/  cos/  *  |  f/3  cot/] 

+  am^f3/ esc/-  4f/cos/ -  £f3/3csc/  cot/  -  ^Isir/cos/ 

♦  f3sin/]  +S5  p'/  +  :Vf3sir/  cos/  +  f/3csc/] 


The  i°  term  of  the  above  series  has  been  maximized  and  this 
value  is  shown  in  Appendix  I,  Part  A..  The  error  in  the  azimuth  °( 
which  would  result  from  the  omission  of  the  f3  term  of  has 

been  recorded  in  Appendix  I,  part  A. 
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In  the  case  of  short  geodetic  lines  (lines  shorter  than  160 
miles)  when  the  values  of  §t  A -L,  etc.  of  the  series  above  become 
small,  it,  is  necessary  to  use  floating  point  notation  in  order  to 
insure  greater  decimal  accuracy. 

The  alternative  formulas  for  sin  cot<^i^j»  and  cot  o^-l*  which 
are  given  below  are  recommended  for  short,  lines.  (They  may  also  be 
used  for  long  lines). 

sin  "^sinL  cos^2^  +  ^sin(f^2-/^l' +  2cos/4  sin/(^sin"  ^  J 

cot<^l_2  -  Lsln(/?2-A)  +  co  sir/^(l  -  cos  A 

cosy^2  s^n  ^ 


cot  <^2-1 


(sin</$2“/^l)  "  c°4  sit  A\  1  -  cos  A)) 


where 


- 11 - - 

cos sir.  / 

2  -  ~  B2-B]_  +  jn(sin2?,;]_-s3r.2B2)”~(sirljli-sinl4E2) 

+  -  (sin6B1-sin6l2)|' 


1  -  Vl-e‘ 


1  +  \[l-ed 

(3^  and  r2  as  previously  defined  are  the  see  -r aside  latitudes  of  pointr 
1  and  ?,  respectively).  t.V.e  our- Irani  cri  t.eri  <  1.  or.  pope  20. 

LlST-lCi  V. 

IT.  ELECTRONIC  COKPnTif;  lUThCO  NCR  SCVj""10a  Of  TlhECT 


'hi' .Tie  r.RCdLfcl.,  -vi?ED  O’. 


h:  i«v’  v 


The  formulas  f.  riven  tel  j”  for  i  hu  -'ol  dp  r,  rt*  the  direct  geodetic 
iroblom  a :'u  in4,  out. 3  i  *  ri.  *rf  ly  for  ol-.chros ;  c  computer  pro>  raming . 
in-wo;,  tic  :  a;  also  \  e  use  I  i>r  0 ut. ■■  ti  nr.  h;  nears  of  desk  cal- 


CU'at CVS. 


TJ’P  computation  form  for  this  method  is  as  follows s 
''riven:  *  pen  rnphic  latitude  and  lonpitude,  respect:  v  ■ 

of  any  point  1. 

®  “  azimuth  clockwise  from  north, and  distance, 
respectively,  to  any  other  point  2, 
keqrired:  Bj,  1<2  and  (Latitudes  and  lonpit^%js  consider 

(+)  north  and  east,  (-)  south  and  west). 


/o  -  [/a] 


♦  ai  ['  f~2  S±n^s  j 

"h  ^  /s  +  sin/s  cos/s  j 

al°  |jr~  sin/s  cos>rsJ 


?  /J-  iae'k  -  rT  rT  e'^ 

+  *1  pp-  /s  “  igp  sin/s  cos&  - 


/s  C0S?^ 


+  ,ir/r 


^3oT  ) 


+  alml^  T —  sin$s  +  Z—  /so°S/s  "  -  Slni^s  cos  i.i  I 

IdlA.  .  =  f"f£s  ]  +  al[~sin/s] 

fco T  L  J  L  ^  J 

+  mx  [Z1  /g  -  Sin/S  cos^Cs  ] 
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where : 


ao>  -  semi -major  and  semi -minor  axes,  respectively,  of 
sr  >her  odd 

f  =  spheroidal  flattening  =  ^1  -  j 


r 


second  eccentricity  squared  =  (a2  -  b^)  a  b2 

O  O  0  G 


number  of  seconds  in  one  radians  =  206 ,26)4, 80625 
j?3  =  (S  «■  bQ)  radians 

tan  $  -  (tan  B)  (1-f)  when |R  Ii5°  or  cot  when 

J  (1  -f) 

|B}> 

cos/(3o  =  cos/^  sinCX^P 

g  ■  cos  A  coso<'1_2 

al  *  (1  +  i_£  sin2/^)  {sin2/$1  cosjf3+  g  sin/^  sin  $,) 

ml  ■  (1  ♦  f?~-  sin2^-,)  (1  -  cos2/3  ) 

U  2 

sin/3  2  ”  sin/^i  cos  §0  +  g  sin  $0 

ccs/o  =  +  >f(cos$oy-  *  {iviXiXl  sin  /0  -  g  cos/0)? 

The  quantity  under  the  radical  and  its  root  must  be  computed  by  ;'V 
^r.  :<c-.-.n  :  !*  ,J- - -.'-f.  _  v;'-  n  '-'••'■pi  -’ulx  nt  1-mvc  chsclu4 

i  A  =  ccs^ 


latitude,,  jf  sin/ 


tan  P2  ~ 


or  co 


,  whichever  has  the 


cos  A?  ' c 

smaller  absolute  value. 

Obtain  tan  (or  cot)  of  B2  from  its  relation  to  tan  (or  cot)  of  A 

1  c.  • 

Obtain  B2>  which  ranges  from  -90°  through  +90°  and  takes  the  sign 
of  its  tan  (or  cot). 

coW2_l  K  (cos cos  jJ0  -  tan/^x  sin  f0)  *  sin 
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When  j  |  >1,  divide  result  into  1  fco  obtain  tan»<2-l  instead. 

(,'r.it  4 here  Ir.st  fee  ?f.noc  for  meridional  arcs) 


Quadrant  of  ofp_  -± 

If  (0%<AU2<180°) 

....  and  cot  (or  tan)  °fo<2-l  is  (+)  or  (-), 
o<  ?  ^  is  in  Quad.  Ill  or  IV,  respectively. 

If  (180^«p<1w2<360o) 

....  and  cot  (or  tan)  of^p,,^  is  (+)  or  (-), 
c/p  i  as  in  Quad*  I  or  II,  respectively. 

For  meridional  arcs,  enter  the  above  table  with  the  sign  of  the 


J  numerator  of  cot  o<2-l>  and  reference  angle  0°* 

cot  A  *  (cot  cos  A  -  cos  °<1  _2  sin  A)  *•  sino<2_2 
■•'hen  |  cot  A  j  >1,  divide  result  into  1  to  obtain  tan  A  instead. 


(>;rit  ■  •  Z.'-r-  -w  13  \  '  for  ;wricio;.nl  arcs) 

Quadrant  and  Sign  of  ^  J 

When  JjL  is  in  Quad. I 
or  II  (l80°  included) 

When  J^o  is  in  Quad. Ill 
or  IV  (180°  excluded) 

and 

(o^c^i.giieo0) 

...  then  if  cot  (or 
tan)  of  A  is  (+)  or  (-) 
A  is  in  Quad.  7  or  II , 

respectively. 

...  then  if  cot  (or 
tan)  of  A  is  (♦)  or  (-) 
A  is  in  Quad.  Ill  or  IV. 
respectively. 

and 

( 180°<  o<1-2<360°) 

...  then  if  cot  (or 
tan) of  A  is  (+;  or  (-) 
the  assoc,  angle  is 
in  Quad.  Ill  or  IV, 
respectively,  and  A 
is  obtained  by  sub~ 
tractin^  360® 

...  then  if  cot  (or 
tan)  of  A  is  (+)  or  (-) 
the  assoc,  ancle  3s 
in  Quad.  I  or  II, 
respectively,  and  A 
is  obtained  by  sub¬ 
tracting  360° 

For  meridional  arcs,  enter  the  above  table  with  the  sign 
of  the  numerator  of  cot  A,  and  reference  angle  0°. 


hp  m  +  L 

If  J  hp  |  >  180°,  modify  L>  by  adding  or  subtracting  360° 
according  to  whether  it  is  initially  negative  or  positive. 


1?.  Extension  of  Tories  of  the  Direct  Geodetic  IToblem  For  Greater 


The  e'°  term  of  the  preceding  series  has  been  derived  and 

numerically  maximized  in  order  to  determine  the  error  in  the  §Q 

series  which  would  result  from  an  omission  of  the  e term.  This 

maximum  value  is  given  in  Aprernix  I,  Tar*-  A. 

Likewise,  the  f I.~ T— ^  )  series  has  been  extended  to  include  the 

Vcos^, 

f3  term.  A  maximum  numerical  value  is  given  in  '  r  p  j ,  .  ... 

The  errors  which  the  omission  of  the  e  ’6  term  of  and  the  f3 
term  of  ^ — —  )  could  finally  produce  in  Bp,  Lg  and  **2-1  are  also 

Pc-oa/Ao 

shown  in  Appendix  2,  Part  A , 

The  e'*  term  of  the  $0  series  is  as  follows: 

|e'5  term  of  »  aj3  |^e'6  sin^cos2/s  «-  sin^J 

+  al2m1|^|  e'%  •  ^  e'6/s  cos 2fs  “  ~e  ,6  sinj^  cosjtf” 


x  1L.6 


■e 1  sm 


>/s  cos3/sJ 


36 


1 


In  the  above  equation  t]  =  sin2^!  and  all  other  quantities 
are  the  same  as  defined  on  page  33. 

As  /«,  approaches  180°,  each  of  the  two  terms  containing 
caqj^  in  the  series  above  anproaches  infinity.  However,  they  may 
he  combined  into  a  single  finite  term  as  follows: 

al  "?  (-^  §£'  cos2/s  cac/s)  ♦  n l  t±  (S^L  /s2  cos^  CBoft)  - 

raiy(S4|  h2  cos/s)  *  where  Y  •  (si n2/$1  sinjQ  -  g  siny^i  cos/s). 


The  series  for 


/cos/? 


extended  through  the  fJ  term  is  as 


follows: 


l^-  M  •  -N  f2  +  2f3)  sin/s| 

+  ml^f2  +  f3)/s  +  (~  “  £3)sin/s  cos/E j 

+  a^2  Uf3  sin^Tg  cosj^sJ 

*  31^1  ^  -|  f3  sin^«  |i'3/s  cos/s 

+i4'3sin^s  cos2^  J 

*  "i2  [  -|f3A  *'zr¥s  »■% 

“  f3sin/scos3/5  +  ^f3sin/s  cos^J 
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